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Abstract: It is well accepted today that diffusive acceleration in shocks results to the cosmic ray spectrum 
formation. This is in principle true for non-relativistic shocks, since there is a detailed theory covering a 
large range of their properties and the resulting power-law spectrum, which is nevertheless not as efficient 
to reach the very high energies observed in the cosmic ray spectrum. On the other hand, the cosmic ray 
maximum energy and the resulting spectra from relativistic shocks, are still under investigation and debate 
concerning their contribution to the features of the cosmic ray spectrum and the measured, or implied, 
cosmic ray radiation from candidate astrophysical sources. Here, we discuss the efficiency of the first 
order Fermi (diffusive) acceleration mechanism up to relativistic shock speeds, presenting Monte Carlo 
simulations. 



Introduction 

The notion that collisionless plasma shocks can 
accelerate charged particles into high energies re- 
sulting into the observed cosmic ray spectrum, has 
been around for many years and by now it is widely 
accepted that this process can account for the ma- 
jority of the energetic, non-thermal particle distri- 
butions that we infer in various astrophysical envi- 
ronments. The work of the late 70s by a number of 
authors (e.g. [5], [1], [2]) established in principle 
the basic mechanism of particle diffusive accelera- 
tion in non-relativistic shocks. Since then, consid- 
erable work has been done analytically and numer- 
ically on the subject, however questions still need 
to be answered on this acceleration mechanism es- 
pecially at relativistic shock speeds. 

In the present paper we will present simula- 
tions from non-relativistic up to relativistic shock 
speeds, in order to test the acceleration efficiency 
of the first order Fermi mechanism. More pre- 
cisely, these simulations include sub-luminal and 
super-luminal shock acceleration with an aim to 
draw conclusions on the primary spectra of the 
relevant astrophysical sources such as Supernovae 
Remnants, Active Galactic Nuclei hot spots and 
Gamma Ray Bursts. 



Numerical simulations 

The purpose of using a Monte Carlo code is to sim- 
ulate the particle transport equation, in our case 
for mild and relativistic shock velocities. The ap- 
propriate time independent Boltzmann equation is 
given by 

tvt/ , \ 9 f d i ' \ 

nv+vp)- = -\ c 

where V is the fluid velocity, v the velocity of the 
particle, V the Lorentz factor of the fluid frame, p 
the cosine of the particle's pitch angle and ^ | c the 
collision operator. The frames of reference used in 
the code are the local fluid frame, the normal shock 
frame and the de Hoffman-Teller frame. Pitch an- 
gle is measured in the local fluid frame, while the 
value x is the distance of the particles to the shock 
front, where the shock is assumed to be placed at 
x = 0. Specifically, in order for the above equa- 
tion to be solved by the Monte Carlo technique, 
we assume that i) the collisions represent scatter- 
ing in pitch angle and ii) the scattering is elas- 
tic in the fluid frame. We begin the simulation 
by injecting 10 5 particles far upstream. First or- 
der Fermi (diffusive) acceleration can be simulated 
provided that the particles' guidance center under- 
goes multiple scatterings with the assumed mag- 
netized media and in each shock crossing the parti- 
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Figure 1 : Two sub-luminal shock spectra, one for 
a non-relativistic shock of a gamma factor equal to 
1.0001 (i.e. y=0.01c) and the other for a gamma 
equal to 5.0 (i.e. V=0.98c) with an inclination an- 
gle of 89 and 39 degrees respectively. One can 
see the almost perfect superposition of the spectra, 
showing the equivalence of acceleration efficiency 
between the non-relativistic perpendicular shock 
and the relativistic sub-luminal oblique one. 

cles gain an amount of energy. The mean free path 
A is calculated in the respective fluid frames (i.e. 
upstream or downstream) assuming a momentum 
dependence to the mean free path of the particle 
such as, A = A Q • p, where A D = 10 • r g , which 
is related to the spatial diffusion coefficient, k. 
Specifically for oblique shocks, in the shock nor- 
mal which lies in the x direction, the diffusion co- 
efficient k, equals to k\\ cos 2 ip, where k\\ — Xv/3 
and ip is the inclination of the shock to the mag- 
netic field lines. A pitch angle diffusion algorithm 
is used (e.g. see [6]). All results are given at the 
downstream side of the shock frame. In general, 
flow into and out of the shock discontinuity is not 
along the shock normal, but a transformation is 
possible into the shock frame to render the flows 
along the normal (e.g. [3]) and for simplicity we 
assume that such a transformation has already been 
made. During our simulations continuous Lorentz 
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Figure 2: Super-luminal shock spectrum for a 
gamma factor of 10 at 89 degrees. 



transformations are performed from the local fluid 
rest frames to the shock frame to check for shock 
crossings. 

Furthermore, for the super-luminal shock condi- 
tions (see [6]) we consider a helical trajectory mo- 
tion of each test-particle of momentum p, in the 
fluid frames upstream or downstream respectively, 
where the velocity coordinates (v x ,v V7 v z ) of the 
particle are calculated in a three dimensional space. 
In principle, we follow the helix trajectory of each 
particle in time t, where t is the time from de- 
tecting the shock intersection at x, y, z. Further 
details of this method can be found in the work 
mentioned above. All particles leave the system if 
they escape far downstream at the spatial bound- 
ary Tb, or alternatively, if they reach a specified 
maximum energy E max , for computational conve- 
nience, even though other physical parameters de- 
scribing particle escape or energy loss would need 
to be taken into account in more realistic situa- 
tions. The downstream spatial boundary required, 
can be initially estimated from the solution of the 
convection-diffusion equation in a non-relativistic, 
large-angle scattering approximation in the down- 
stream plasma which gives the chance of return to 
the shock as, exp(— V^Tb jxi)- In fact, we have per- 
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Figure 3: Two super-luminal shock spectra for a 
gamma factor of 500, one for an inclination an- 
gle of 50 degrees and the other for an inclination 
of 89 degrees. The spectra reach a maximum of 
10 14 eV for protons. This cut-off in the energy 
is due to the stream flow, sweeping the particles 
downstream limiting their chance of returning back 
to the shock. 



formed many runs with different spatial boundaries 
to investigate the effect of the size of the accel- 
eration region on the spectrum, as well as to find 
a region where the spectrum is size independent. 
We note that in the pitch angle diffusion case, the 
inherent anisotropy due to the high downstream 
sweeping effect may greatly modify this analytical 
estimate. 

Results 

In this section we present results of the simulations 
shown in the figures 1-4. 

Two sub-luminal shock spectra are shown in fig- 
ure 1, one for a non-relativistic shock of a gamma 
factor equal to 1.0001 (i.e. V=0.01c) and the other 
for a gamma equal to 5.0 (i.e. V=0.98c) with an 
inclination angle of 89 and 39 degrees respec- 
tively. One can see the almost perfect super- 



Figure 4: A super-luminal shock spectrum for a 
gamma factor of 900 at 75 degrees. One sees 
(simulated spectra for other shock inclinations are 
comparable) the smooth power-law formation of 
the spectrum. Nevertheless, even for such high 
Lorentz factor the super-luminal shocks do not 
seem to be as efficient as the sub-luminal ones (e.g. 
[6]), since the simulations showed that almost all 
of the particles are advected downstream just after 
a shock cycle. 



position of the spectra, showing the equivalence 
of the acceleration efficiency between the non- 
relativistic perpendicular shock and the relativis- 
tic sub-luminal oblique one. The above correlation 
of efficiency between non-relativistic perpendicu- 
lar shocks and mild relativistic oblique ones, ap- 
plies to many other shock spectra we simulated, as 
long as the non-relativistic shocks are nearly per- 
pendicular. For the case shown here, the maxi- 
mum energy corresponds to ~ 10 7 5 GeV for pro- 
tons. Adding to the above, [9] have shown that 
non-relativistic perpendicular shocks, are efficient 
accelerators for cosmic rays reaching energies as 
high as ~ 10 17 eV, since the acceleration in those 
shocks seems to be faster under certain diffusion 
conditions [4]. Also as [6] have shown, sub- 
luminal highly relativistic shocks are the most ef- 
ficient accelerators comparing to the above, result- 
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ing in an average cosmic ray energy gain of <~ T 3A 
after two complete shock cycles. For this case, T is 
the boost factor of the shock, while in the simula- 
tions the test-particles are considered already rela- 
tivistic with an initial Lorentz gamma, 7 = T+100. 
One further result of the simulations, is that for 
oblique shocks below a Lorentz factor of T=5 (i.e. 
V=0.98c) as higher the inclination of the shock as 
greater its efficiency to accelerate particles up to 
10 17 eV (i.e. protons, iron nuclei). 
In figure 2 we show a super-luminal shock spec- 
trum for a gamma factor of 10 at 89 degrees. In 
figure 3 one sees two super-luminal shock spectra 
for a gamma factor of 500, one for an inclination 
angle of 50 degrees and the other for an inclination 
of 89 degrees. Simulations show that the inclina- 
tion of the shock in the super-luminal cases does 
not change the results. The cut-off is due to the 
flow, sweeping the particles downstream limiting 
their chance to return back to the shock, diminish- 
ing the chances for these shocks to be efficient cos- 
mic ray accelerators. In figure 4, a super-luminal 
shock spectrum for shock gamma factor of 900 at 
75 degrees is shown. One sees the smooth power- 
law formation of the above spectra (in contrast to 
the sub-luminal spectra formations in [6]). Never- 
theless, even for as high Lorentz gamma as 900, 
the accelerated particles do not reach very high en- 
ergies (e.g. 10 14 eV for protons) since almost all 
of them are advected downstream, after their helix 
trajectory performs a shock-cycle crossing. All the 
simulated super-luminal shock spectra can be well 
fitted by a power-law with a spectral index rang- 
ing between ~ 2.0 — 2.3 in contrast to the work of 
[6] concerning sub-luminal relativistic shocks. The 
latter spectra appear flatter, at the highest speeds 
with a characteristic plateau-like structure, since 
first, a highly relativistic shock catches the parti- 
cles up within less a shock cycle, which does not 
allow sufficient time for isotropisation and second, 
because relativistic shocks are sensitive to the ap- 
plied particle scatter model [8] . 

Conclusions 

We discussed the mechanism of diffusive shock ac- 
celeration, presenting simulations in order to test 
the efficiency of sub-luminal and super-luminal 
relativistic shocks in a quantitative comparison to 



non-relativistic ones. It is certain that super-luninal 
shocks are not as efficient first order Fermi cosmic 
ray accelerators as the sub-luminal ones. Never- 
theless, the super-luminal shocks can be well fit- 
ted by a power-law with a spectral index rang- 
ing between ~ 2.0 — 2.3 in contrast to the rela- 
tivistic sub-luminal ones which give flatter spec- 
tra. Specifically, the super-luminal shocks are ef- 
ficient in accelerating particles up to a maximum 
energy of around 10 14 eV (i.e. protons or iron nu- 
clei). On the other hand, the sub-luminal relativis- 
tic oblique shocks are better accelerators concern- 
ing a maximum energy above 10 18 eV, as shown in 
similar works as well, but the spectral features are 
not as smooth as in super-luminal ones (which are 
given by a 'clean' power-law), indicating a strong 
connection to the kinematical details of the parti- 
cle diffusion while crossing the relativistic shock 
front. Interestingly, it is shown that perpendicu- 
lar non-relativistic shocks seem to be as efficient 
accelerators as sub-luminal mild relativistic ones. 
There is work under way concerning the initial pa- 
rameters which could affect the shock behaviour in 
highly relativistic speeds. 
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